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1.  Introduction 

Today,  second-  and  third -generation  wave  models  are  used  for  many  applications 
around  the  world.  From  the  design  of  coastal  structures  to  proper  regional  sediment 
management  and  from  estimates  of  mixing  in  the  upper  ocean  to  the  efficient  the  routing 
of  ships  across  oceans,  the  accuracy  of  these  models  has  critical  importance.  However, 
following  the  development  of  the  WAM  model  in  the  mid  to  late  1980’ s,  there  has  been 
little  effort  focused  directly  on  wave  model  development.  Instead,  most  research  in  this 
area  has  involved  re -calibrating  various  source  terms  to  obtain  better  fits  to  observations 
(satellite  and  in  situ).  Two  likely  reasons  for  the  lack  of  focus  in  this  area  are  1)  the 
general  acceptance  outside  of  the  wave  community  that  the  WAM  model  physics 
provided  a  good  “detailed -balance”  description  of  the  wave  generation  process  and  2) 
wave  researchers  cannot  agree  on  necessary  modifications  to  the  WAM  physics. 

The  purpose  of  this  paper  is  to  revisit  the  general  framework  of  the  physics  of  wave 
generation  and  attempt  to  formulate  a  concept  for  wave  generation  that  seems  to  “fit” 
some  important  constraints  implicit  in  observations  over  a  wide  range  of  generation 
scales  (time  and  space).  Toward  this  end,  data  from  three  very  different  sites  will  be 
used: 

1.  Lake  George  -  a  small,  shallow  site  located  in  Southeast  Australia  (Figure  1); 

2.  FRF  Gage  #630  -  a  coastal  site  located  in  the  Atlantic  Ocean  in  a  depth  of 
about  18  meters  of  water  approximately  5  km  off  the  coast  of  Duck,  North 
Carolina  (Figure  2);  and 

3.  NDBC  Buoy  46035  -  an  open -ocean  site  located  in  deep-water  in  the  Bering 
Sea  (Figure  3). 

These  sites  cover  a  wide  range  of  practical  wave  generation  scales,  from  shallow  to  deep 
and  from  very  small  fetch  to  extremely  large  fetch. 


2.  Theoretical  Framework 

Ever  since  the  mid  1970’s,  most  simulations  of  the  evolution  of  wave  spectra 
under  the  action  of  the  wind  have  included  the  effects  of  wave-wave  interactions. 
Unfortunately,  this  nonlinearity  effectively  couples  the  entire  wave  spectrum  into  a  single 
dynamic  system,  making  it  diffic  ult  to  assess  the  precise  roles  of  different  forcing 
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mechanisms.  In  an  attempt  to  avoid  this  problem,  we  postulate  a  theoretical  framework 
for  wind -wave  spectra  during  active  wave  generation  as  shown  in  Figure  4.  In  this 
Figure,  four  spectral  regions  are  defined,  based  roughly  on  the  hypothesized  controlling 
forcing  mechanisms  for  each  region. 

Region  1,  the  spectral  peak  region,  is  bounded  at  its  upper  limit  by  frequency  fo, 
which  is  a  frequency  at  which  there  is  no  net  nonlinear  flux  of  energy  and  extends  to  zero 
at  its  lower  limit.  Since  there  is  no  net  flux  of  energy  into  or  out  of  Region  1  due  to 
wave-wave  interactions,  the  net  gain  or  loss  of  energy  in  this  region  depends  only  on 
external  sources  (wind  input  and  wave  breaking).  The  role  of  wave- wave  interactions 
region  will  likely  relate  to  down-shifting  of  the  spectral  peak,  due  to  the  strong 
asymmetry  of  this  source  term  within  this  region. 

Region  2,  the  transition  region  from  Region  1  to  the  equilibrium  range,  is 
bounded  by  fo  at  its  lower  limit  and  by  fgq  (the  low-frequency  limit  of  the  equilibrium 
range)  at  its  upper  boundary.  This  region  represents  an  area  in  which  nonlinear  energy 
fluxes  become  less  and  less  influenced  by  the  proximity  to  the  spectral  peak  (  and  low- 
frequency  cut-ff)  as  feq  is  approached.  Since  there  is  zero  net  energy  flux  across  fo,  we 
will  hypothesize  here  that  the  primary  mechanisms  controlling  the  shape  of  this  region 
are  net  wind  input  and  wave-wave  interactions.  Under  the  assumption  that  wave 
breaking  is  very  small  in  this  region,  a  positive  flux  of  energy  (i.e.  a  flux  directed  toward 
high  frequency)  must  exist  at  Ij-q,  in  order  to  compensate  for  the  energy  gain  in  this 
region. 


Region  3,  the  equilibrium  range,  has  been  the  subject  of  many  studies  (Zakharov 
and  Filenenko,  1966;  Kitaigorodskii,  1983;  Resio,  1987;  Resio  et  al,  2001).  It  is  assumed 
here  that  this  region  extends  from  fgq  at  its  lower  limit  to  some  high  frequency  at  which 
dissipation  begins  to  play  a  major  role  in  the  spectral  energy  balance.  In  the  initial 
derivation  by  Zakharov  and  Filenenko  (1966),  it  was  assumed  that  an  (where  (O  is 
radial  frequency)  slope  would  exist  in  this  region  ( i.e.  E(0i))~t0"4,  where  E((D)  is  energy 
density)  in  the  absence  of  any  external  sources  and  sinks.  However,  since  we  know  that 
wind  input  exists  at  the  spectral  peak  and  at  higher  frequencies,  it  is  extremely  unlikely 
that  wind  input  is  identically  zero  across  the  entire  equilibrium  range.  Consequently,  it 
would  appear  that  the  requirement  of  no  input  should  be  relaxed  to  a  sufficiently  small 
input  such  that  the  equilibrium  slope  is  not  moved  far  from  the  (0'"^  slope  observed  in 
most  deep-water  data  sets.  Recently,  Resio  et  al  (2001)  showed  that  a  more  general  form 
for  the  equilibrium  range  could  be  written  in  term  of  a  wavenumber  spectrum, 

F(k)  =  ^u,g-‘''-k-='^  (1) 

where  E(k)  is  the  energy  density  in  wavenumber  space,  a4  is  a  universal  constant,  Ug  is  a 
scaling  velocity  ,  g  is  gravity,  and  k  is  wavenumber. 


We  assume  that  all  three  traditional  deep-water  source  terms,  wind  input  (Sin), 
nonlinear  wave-wave  interactions  (Sni),  and  wave  breaking  (Sds),  occur  in  all  regions  of 


the  spectrum  to  some  degree.  For  a  stable  equilibrium  to  exist,  the  sum  of  all  three 
sources  must  equal  zero.  In  the  equilibrium  range,  if  the  sum  of  the  wind  and  wave 
breaking  source  terms  are  not  zero,  Sni,  which  is  equal  to  the  divergence  in  the  nonlinear 
energy  fluxes,  must  be  non- zero  to  maintain  a  net  zero  balance.  In  other  words,  we 
would  have 


^  =  +  (2) 

dCO 

where  Fe  is  the  flux  of  energy  through  the  spectrum.  If  wind  input  and  wave  breaking 
exactly  cancel  each  other  in  this  range,  equation  2  would  still  hold.  It  would  only  mean 
that  the  divergence  would  be  identically  zero  over  this  range  and  the  fluxes  should  be 
exactly  constant. 

Nonlinear  energy  fluxes  transfer  energy  from  low  to  high  frequencies  and  from 
high  to  low  frequencies  simultaneously.  In  some  previous  papers,  the  term  flux  has  been 
used  to  describe  the  positive-directed  fluxes  (fluxes  from  low  to  high  frequency),  while 
the  negative-directed  fluxes  (fluxes  from  high  to  low  frequency)  are  termed  inverse 
fluxes;  however,  this  terminology  will  not  be  adopted  here.  We  shall  take  the  meaning  of 
Fe  to  be  the  net  difference  of  the  fluxes  in  both  directions,  i.e. 


r^  =  r,"+r;  0) 

where  the  superscripts  “-r”  and  denote  fluxes  toward  higher  and  lower  frequencies, 
respectively.  According  to  essentially  all  recent  theoretical  treatments  of  Sni,  nonlinear 
fluxes  within  the  equilibrium  range  tend  to  force  the  spectrum  toward  an  co  '*  spectral 
form  in  the  absence  of  other  significant  source  terms.  However,  as  will  be  shown  here, 
even  for  typical  wind  inputs,  the  deviation  from  an  (o  '*  spectral  form  may  be  quite  small, 
Within  the  equilibrium  range,  the  flux  of  energy  from  lower  to  higher  frequencies  is 
given  by  (Resio  et  al.,  2001)  in  a  form  equivalent  to 

r;=c.,7i'p-'  (4) 


where  Cni  is  a  dimensionless  coefficient  depending  weakly  on  angular  spreading  and 
proximity  to  the  spectral  peak  and  (3  is  a  wave  steepness  parameter  equivalent  to 
a4u/(2g  )  in  equation  2.  To  a  first  approximation,  the  divergence  of  the  flux  can  be 
written  as 


dco  dco 


(5) 


where  ^  is  a  compensated  spectral  density  of  the  form 

a^u 


(6) 


From  these  arguments  it  ean  be  seen  that  a  20%  inerease  in  speetral  densities  in  the 
equilibrium  range  would  result  in  over  a  70%  increase  in  the  energy  flux  rates  due  to 
nonlinear  interactions.  Thus,  the  nonlinear  interactions  should  force  a  very  strong 
tendency  toward  the  (0'"^  spectral  form  even  if  wind  input  is  occurring  in  this  region  of  the 
spectrum. 

Region  4,  the  dissipative  range,  is  the  portion  of  the  spectrum  in  which  energy, 
fluxed  toward  high  frequencies  via  nonlinear  interactions,  is  lost  due  to  either  turbulence 
(breaking)  or  viscous  effects.  The  low-frequency  limit  to  this  range  should  be 
coincident  with  the  point  where  the  equilibrium  range  slope  shifts  to  a  steeper  slope  (O)'" 
where  n>4).  This  may  occur  over  a  range  of  frequencies  and  different  researchers  have 
postulated  different  points  at  which  this  is  expected  to  occur.  For  the  purpose  of  this 
paper,  we  will  just  assume  that  a  high-frequency  energy  sink  exists,  one  that  is  capable  of 
removing  all  energy  fluxed  into  this  region  of  the  spectrum.  Since  this  spectral  region  is 
located  at  relatively  high  frequency,  the  exact  shape  of  the  spectrum  in  this  region  does 
not  affect  most  practical  applications. 


3.  Data 


Data  selected  for  analysis  here  include  one  deep  ocean  buoy  maintained  by  the 
National  Data  Buoy  Center  (NDBC,  htty ://www. ndbc. noaa.gov).  one  shallow  ocean 
sensor  maintained  by  the  U.S.  Army  Engineer  Research  and  Development  Center’s  Field 
Research  Facility  (FRF,  http://www.  frf.  usace.  army.mil),  and  one  site  on  the  east  side  of 
Lake  George,  Australia.  Anemometers  were  located  directly  above  wave  sensors  at  all 
sites  except  for  FRF  #630,  where  winds  were  measured  at  the  landward  end  of  the  FRF 
pier. 


Four  dimensionless  parameters  can  be  considered  relative  to  site  selection: 
relative  depth  of  the  waves  in  the  equilibrium  range,  ,  relative  wave  height  /  h  , 

wave  steepness  ,  and  inverse  wave  age  u^/c^  (often  characterized  as  u^QlCp). 

The  first  of  these  has  an  approximate  theoretical  lower  limit  ~  0.7  below  which  the 

wavenumber  dependence  of  the  equilibrium  range  may  change  form,  due  to  depth -related 
changes  in  the  behavior  of  coupling  coefficient  (Resio  1987).  Zakharov  (1999)  has 
shown  that  for  k^^h  <  0.3  ,  the  spectrum  should  asymptotically  tend  toward  a  k~‘^'^  form. 

The  subscript  “eq”  has  been  added  to  his  original  “kh”  parameter  to  reflect  the  fact  that 
his  theoretical  argument  is  valid  only  within  the  equilibrium  range  and  not  in  the  spectral 
peak  region.  With  parameters  defined  in  Sections  3.2  and  3.3,  data  from  the  above  sites 
have  relative  depths  in  the  range  0.7  <  k^Ji, ,  which  pushes  the  lower  limit  somewhat. 

The  data  extends  into  very  deep  water  at  most  deep -ocean  sites. 


4.  Analysis 


Figure  5  shows  a  plot  of  the  band- averaged  valued  for  ^  as  a  function  of  f/fp,  for 
all  three  data  sets  described  later  in  this  paper.  In  this  Figure,  the  subscript  “p”  denotes  a 
value  at  the  spectral  peak.  Also  shown  in  this  Figure  are  the  curves  for  equivalent 
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in  deep  water)  and  k’  ((O'  equivalent  in  deep  water)  spectra,  which  clearly  do  not  fit  the 
data  very  well.  Even  though  some  random  deviations  and  some  possibly  periodic 
variations  in  the  value  of  ^  are  seen,  the  general  tendency  is  for  ^  to  remain 
approximately  constant  in  aU  of  the  data  sets. 

Since  we  are  interested  here  in  an  integral  balance  of  energy,  rather  than  a 
detailed  frequency-direction  balance,  we  shall  not  attempt  to  build  upon  different  source 
term  forms;  rather,  we  shall  examine  six  types  of  wind-input  scaling.  These  can  be 
written  in  the  form: 


(('s, „(/)#-«;  (7) 

and 

D's, „(/)#-«;<■,  (8) 

where  Cp  is  the  phase  speed  of  waves  at  the  spectral  peak,  Ua  is  a  scaling  parameter  with 
units  of  velocity,  taken  here  to  be  1)  the  wind  speed  at  a  reference  level  of  10  meters  (uio) 
;  2)  the  friction  velocity  (u*  =x^^^/Pa,  where  x  is  the  wind  stress  and  Pa  is  air  density);  and 
3)  the  wind  speed  at  a  fractional  height  of  the  wavelength  of  the  spectral  peak  (u;i,  as 
defined  in  Resio,  et  al.,  1999).  Thus,  the  six  scaling  combinations  investigated  here 
consist  of  three  different  scaling  parameters  for  wind  speed  inserted  into  the  two  different 
scaling  forms  shown  in  equations  7  and  8. 


4.  Results 


Figures  6a-6c  show  a  plot  of  (3  versus 


,  where  Ua,  in  turn,  represents  the  three 


scaling  options  based  on  wind  speed  parameters  only.  Clearly,  this  scaling  option  does 
not  provide  a  unified  scaling  form  for  the  different  data  sets.  Instead,  the  slope  of  the 
small-basin  data  (Lake  George)  is  considerably  shallower  than  the  slope  of  the  open 
ocean  data  sets  (FRF  #630  and  NDBC  46035).  Figures  7a-7c  show  a  plot  of  (3  versus 

,  where  Ua,  in  turn,  represents  the  three  scaling  options  based  on  wind  speed  and 

phase  velocity  parameters  in  the  manner  described  above.  This  scaling  velocity  provides 
a  basis  for  the  data  in  the  slopes  of  each  data  set  are  all  consistent.  Thus,  this  scaling, 
rather  than  a  scaling  based  on  wind  speed  alone,  provides  a  better,  scale -independent 
representation  for  energy  densities  in  the  equilibrium  range. 


Of  the  three  different  velocity  scales  examined,  the  velocity  referenced  at  a  level 
of  0.065  times  the  peak  wavelength  above  the  mean  surface  consistently  provided  the 
scaling  relationship  with  the  lowest  scatter  around  the  regression  line.  In  fact  as  shown  in 
Resio  et  al  (2002),  this  scatter  is  reduced  by  a  factor  of  three  to  four  in  high-quality, 
small-basin  data  sets. 

4.  Discussion 

As  shown  in  Resio  et  al  (2002),  the  findings  here  are  consistent  with  a  constant 
flux  of  momentum  from  the  atmosphere  into  the  wave  field,  which  agrees  with  the  work 
of  Hasselmann  et  al  (1973)  and  Resio  and  Perrie  (1989).  In  fact,  in  the  work  by  Resio 

and  Perrie  (1989)  was  the  initial  basis  for  attempting  the  ulcp  form  as  the  scaling 

parameter  for  (3,  after  we  had  found  that  none  of  the  forms  based  on  wind  speed  along 
gave  a  very  good  result. 

The  best- fit  relationship  for  the  entire  data  set  was  clearly  achieved  with  the 

scaling  parameter.  If  we  return  the  work  of  Resio  and  Perrie  (1989),  we  see  that  they 
included  an  analysis  of  laboratory  data  in  their  scaling  analysis  and  showed  that  a 

MjoCp  form  provided  a  better  fit  to  a  wide  range  of  data  than  did  forms  based  on  wind 

speed  alone;  however,  they  did  not  have  the  extensive  data  shown  here  supporting  that 
conclusion  in  their  paper.  If  we  return  to  the  data  sets  shown  in  their  paper  we  see  that 
the  use  of  u,.  improves  the  consistency  of  the  (3-Ua  relationship  for  lab  data  included. 

These  findings  could  potentially  have  significance  to  at  least  three  several 
practical  applications.  First,  the  interpretation  of  scatterometer  data  is  roughly  based  on  a 
linear  relationship  between  wind  speed  and  measured  return  from  wave  fields.  Since  the 
calibration  of  this  relationship  is  based  primarily  on  open -ocean  data,  it  is  possible  that 
the  actual  relationship  is  more  of  the  form  shown  here.  Second,  wave  models  that  depend 
on  detailed -balance  arguments,  typically  have  been  calibrated  to  reflect  either  an  O)'^ 
spectral  form  or  a  linear  dependence  of  the  energy  levels  in  the  equilibrium  range  on 
wind  speed.  It  is  unlikely  that  either  of  these  detailed  -balance  forms  can  work  well  over 
a  wide  range  of  dynamic  scales.  And  third,  arguments  of  wave  induced  drag  coefficients, 
which  are  strongly  dependent  on  wave  age,  are  not  supported  by  our  findings. 

Another  point  worth  discussing  in  this  paper  is  the  fact  that  considerable  evidence 
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supporting  the  k"  equilibrium  range,  even  into  kh  values  as  low  as  0.6,  is  beginning  to 
accumulate.  This  is  consistent  with  the  expected  behavior  of  spectral  energy  densities  in 
the  equilibrium  range,  based  on  nonlinear  energy  fluxes;  but  it  is  difficult  to  see  how  any 
wave  breaking  mechanism,  which  may  be  coerced  to  depend  on  f"^in  deep  water,  can  be 
made  to  scale  appropriately  to  reproduce  this  result.  This  appears  to  support  arguments 
that  the  primary  energy  losses  within  a  wave  spectrum  are  located  in  at  high  frequencies 
and  not  in  Regions  Till  of  the  spectrum.  It  is  also  consistent  with  physical  arguments 
based  on  the  depth  of  the  “turbulence”  (white-capping)  region  at  the  water  surface  during 


wave  generation,  since  —  « 1 ,  where  5z  is  the  depth  of  the  white-capping  region  and 
Lp  is  the  wavelength  of  the  spectral  peak. 


5.  Conclusions 

The  two  primary  conclusions  based  on  the  analyses  and  data  sets  presented  here 
are  as  follows: 

1.  Energy  densities  in  the  equilibrium  range  depend  on  a  combination  of  wind 
speed  and  spectral  peak  phase  speeds,  consistent  with  the  concept  that  a  constant 
proportion  of  the  momentum  transfer  from  the  atmosphere  enters  the  wave  field. 

2.  Wind  inputs  into  waves  appear  to  scale  better  with  wind  speeds  measured  at  a 
peak  wavelength  scaled  level  above  the  surface  than  with  either  friction  velocity  or  winds 
at  a  fixed  height. 
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Figure  1.  Focation  of  Fake  George  data  site. 
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Figure  2.  Location  of  FRF  #630  site. 
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Figure  3.  Location  of  NDBC  #46035  site. 


Figure  4.  Schematic  of  diagram  for  spectral  energy-flux  regions. 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

Lake  George 


to  6  Hz/^, 


I  I  L  I  I  I  I  I 


"I" 'T' 


I'  '  'f"v'  r" 


•'•"T  •' 


I]  '  ^ 

. . 


FRF  #630 


NDBC  46035 


“’-fifiiiiiH'lift}}  H  fi 


Q  I - 1 - 1 - 


ijl  Jl  11 

I  ...  I 


Figure  5.  Behavior  of  spectral  density  at  data  site. 


Figures  6a  -  c. 


Figures  7a  -  e. 


